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Abstract. Recent observations revealed that small magnetic elements abundant at the solar 
surface move poleward with a velocity which seems to be lower than the plasma velocity U. 
Guerrero et al. (2011) explained this discrepancy as a consequence of diffusive spreading of the 
magnetic elements due to a positive radial gradient of \Ue\. As the gradient's sign (inferred 
by local helioseismology) is still unclear, cases with a negative gradient are studied in this 
paper. Under this condition, the velocity of the magnetic tracers turns out to be larger than the 
plasma velocity, in disagreement with the observations. Alternative mechanisms for explaining 
^0 them independently are proposed. For the turbulent magnetic pumping it is shown that it has 
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to be unrealistically strong to reconcile the model with the observations. 
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d 1. Introduction 

[ Amongst the axisymmetric constituents of the plasma flow in the Sun's convection 
^-H ' zone, the meridional circulation (MC), although being much slower than the differential 
^ ' rotation, is gaining growing interest which arises from its potential importance for the 
solar dynamo. Certainly any advanced non-linear mean-field model of this dynamo will 
have to include the MC to cover the full interaction between mean field and motion. In 
^f-^ particular however, fiux-transport dynamo models depend already on the kinematic level 
• I crucially on the MC as it is the ingredient which allows to close the dynamo cycle. At 
I . the same time it explains naturally both the equatorward migration of active regions and 

the poleward drift of weak poloidal field elements during the solar cycle. 
I . The MC in the Sun is accessible to a variety of techniques out of which Doppler mea- 
^ I surements and helioseismological time-distance or ring diagram analyses tend to agree in 
• . the surface speed profiles with a peak of w 20ni/s at a latitude « 35°. However, attempts 
K% ' to measuring the MC velocity using small magnetic structures as tracers of the flow 
^ . (Komm et al., 1993, Hathaway & Rightmire, 2010) reveal systematic differences from 
" " ' these results. Compared to Doppler measurements (Ulrich 2010) magnetic feature track- 
ing (MFT) speeds can be as much as « 30% smaller. This indicates that the magnetic 
field is not completely frozen in the plasma, which is reasonable given the assumed high 
turbulent magnetic diffusivity 7]^. Due to diffusive spreading of the field, its advection is 
affected by the depth dependence of the MC velocity. On top of advection and diffusion, 
turbulent motions may be contributing by an extra advection term due to turbulent 
magnetic pumping or by the more complicated interaction with dynamo waves. 
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Figure 1. a) Radial profile of [/e(r, 57°) for different values of n. Dashed line: radial pro- 
file of 79, see eq. (|3.1|1 . b) Corresponding tracking velocities, V, for three cases: 9r|(7e| > 
(yellow/continuous) and 9r|(79| < with (dashed) and without (red/dot-dashed) pumping. 
rjT — 10^'^cm^/s throughout. 

In a previous work (Guerrero et al., 2011, hereafter GRBD) we studied numerically 
the kinematic evolution of small bipolar magnetic elements at different latitudes in the 
northern hemisphere. For simplicity we considered turbulent diffusion and advection due 
to a predefined MC as the only magnetic transport mechanisms. In a ID (i.e., surface) 
version of the model, the northern spot acquires an effective velocity which is higher than 
the flow speed, but the southern one travels slower than the fluid. However, on average, 
or considering the center of the bipolar region, its speed clearly fits the fluid speed. This 
result is independent of the value of rjT- In the 2D model, on the other hand, we found 
that the MFT speed differs from the fluid velocity provided that the frozcn-in condition is 
not satisfied. Moreover, we demonstrated that the relevant threshold value of ryx depends 
on the radial gradient of the latitudinal velocity, dr\Ue\, assumed positive. The difference 
between MFT and flow speeds increases with dr\Ue\- Thus, this simple model is able to 
explain the discrepancy between the Doppler and MFT observations. 

Modern local helioseismological techniques allow also to infer the radial variation of 
the MC speed. Unfortunately, time-distance and ring-diagram analyses have not given 
consistent results. The first method yields that at the first few megameters beneath the 
photosphere the latitudinal velocity decreases inwards, i.e., a positive 9r|C/e| (Zhao et al. 
2011). By contrast, the results from the second method indicate a slight inward increase, 

1. e., a negative gradient albeit restricted to middle latitudes. (Gonzales Hernandez et al. 
2006). Such velocity profiles were not considered in GRBD and the aim of this paper is 
to fill this gap as well as to discuss other effects that may infiuence MFT speeds. 

2. MFT speed with a negative radial gradient of \Ue\ 

The method followed here to simulate the evolution of small bipolar magnetic struc- 
tures is described in detail in GRBD. We also use the same notation. Briefiy, we solve the 
induction equation in axisymmetry for the azimuthal component of the vector potential, 
in the variables r and 9. The prescribed meridional flow fed into the model is given 
by eqs. (4-6) of GRBD and has a single cell per hemisphere. In the present paper ijt 
is assumed to be constant across the domain. Bipolar regions (with spot separation 3° 
and depth 0.04i?) at 20 different latitudes from the pole to equator are evolved over a 
time interval of two weeks. Their velocity is computed as an average over this time span 
through the differences in final and initial latitudes (see eq. (8) of GRBD). To obtain 
a negative radial gradient of \Ue\ we preserve its general mathematical form (eq. (5) of 
GRBD), but choose the parameters of this ansatz as fc = 6 and n = 1.4 which results in 
a negative dr\Ue\ within a shallow surface- near layer of roughly 2% of the thickness of 
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Figure 2. Magnetic field lines and r component (color coded) after two weeks propaga- 
tion, rjT ~ lO^^cm^/s. Velocity profile as in Fig. [TJ red curve, a) without, b) with pumping, 
709 = lOm/s, see eQ. (|3.ip 

the convection zone (see Fig. [1^). It is worth mentioning that the stress- free boundary 
condition cannot be satisfied with a negative gradient. With this Ug profile the velocities 
of the magnetic tracers turn out to be larger than the fiow velocity (red/dot-dashed line 
in Fig. [T|d). This result does not come unexpected since the field lines underneath the 
surface are pushed northwards by a flow faster than the surface one. By magnetic tension 
the field lines at the surface are then also accelerated. A snapshot of the magnetic field 
after two weeks of evolution is shown in Fig. [2^. 



3. Additional transport mechanisms 

Given that a negative dr\Ue \ does not allow to reproduce the observations of decelerated 
magnetic tracer motion at the surface, one could be tempted to argue that the results 
of the ring diagram analysis could be incorrect. This conclusion, however, would be 
premature as our model is by no means complete. As mentioned in GRBD, studying 
axisymmetric, i.e., averaged fields requires the inclusion of the full mean electromotive 
force, from which only the turbulent diffusion term has so far been considered. 

What about other constituents? Under anisotropic conditions, turbulence acts like a 
mean advective motion, called turbulent magnetic pumping and is described by a vector 
7. In the northern hemisphere, its 0-component is found to be positive, i.e., equatorwards 
whereas its radial component is mainly negative (downwards) , see |Kapyla et al. 2006a[ 
Being helical, convective turbulence in the rotating Sun exhibits also an a-efFect which 
together with differential rotation in the uppermost 35Mm of the convection zone (the so 
called near-surface shear layer) may well enable a turbulent dynamo of the afl type. As 
dri^ < there and a is expected to be positive in the northern hemisphere, the dynamo 
wave should, according to the Parker- Yoshimura rule, propagate equatorwards, contrary 
to the MC. This scenario is difficult to study as it requires establishing a full nonlinear 
mean-field dynamo model of the surface shear layer. Turbulent pumping, in contrast, can 
easily be introduced as an additive contribution to the MC. We do so employing profiles 



for 7 adapted from Guerrero & de Gouveia Dal Pino (2008) (hereafter GDP08) 
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Similar to GDP08, we set 7or and 706/ such that the maximum of 7e is 2.5 times larger 
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than that of 7r- Q The radial profile of 79, which is the more relevant component for 
this study, is shown with a broken line in Fig. [1^. From direct numerical simulations 
(DNS) |Kapyla et al. 2006a| have estimated that the maximum 7^ should be « 2.5m/s. 
According to eq. p.ip its value at the surface would then be l.llm/s. Given that the 
surface amplitude of jt/gl w 25m/s, this small value cannot significantly modify the results 
found in the previous section. Hence we have increased 700 progressively and found that 
an MFT speed profile that roughly coincides with that of the fiow speed is obtained with 
7oe — 5m/s (i.e., 2.2m/s at the surface or k, 10% of the flow speed). For larger values 
of 7oe the MFT speed is smaller than the fiow speed. In Fig. we show the result 
obtained with the extraordinarily high value 709 = lOm/s (see broken lines). Note that 
the latitudinal profile of 79 (eq. p. II) ) peaks at 63° while the flow profile does at 57°. 
This is interesting because the resulting profile of the MFT speed peaks at a different 
latitude than the fiow, just as it is obtained in observations (see fig. 10 of Ulrich, 2010). 
Fig. [2] compares the morphology of the magnetic field from the simulation with pumping 
(b) to that without (a). Note that the magnetic field is transported further to the north 
in the case without magnetic pumping 

Discussion 

Negative gradients dr\Ug\ at the solar surface, as suggested by helioseismological ring- 
diagram analyses, result in MFT speeds higher than the plasma flow speed. By enhancing 
our flux transport model with turbulent magnetic pumping, we tried to reconcile the 
helioseismological results with the actually oberserved lower MFT speeds. However, high 
amplitudes of the pumping velocity not supported by DNS are needed. We have to 
conclude that either: (i) the negative gradient of |?7e|, suggested by ring-diagram analysis, 
if real, should be smaller than considered here, (ii) the pumping velocities from DNS are 
unrealistically low, which cannot be excluded because the setup in |Kapyla et al. 2006a 
differs from the solar conditions, or (iii) the inclusion of the a effect in the flux transport 
model, enabling a dynamo process in the Sun's surface shear layer, is crucial. For the 
latter, we have to think of a short-wavelength dynamo wave traveling equatorwards. 
Magnetic quenching would cause a corresponding 9 modulation of a, 7 and ryx whereas 
the mean Lorentz force would modulate Ug and il. Here, quenching of 77T is perhaps most 
promising as the value of t^t "decides" to what degree the field is frozen in the fluid. 
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